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Fermionizing the charge sector and bosonizing the spin part in the SU(2) slave-boson theory, we 
derive an effective field theory for dynamics of doped holes in the antiferromagnetically correlated 
spin background, where spin fluctuations are described by an SO(5) Wess-Zumino-Witten (WZW) 
theory while dynamics of doped holes is characterized by QED3 with a chemical potential term. An 
important feature of our effective field theory is the coupling term between valance bond fluctua- 
tions and doped holes. Considering that valance bond fluctuations are deeply related with monopole 
excitations of staggered U(l) gauge fields in the bosonic field theory for spin fluctuations, we demon- 
strate that hole dynamics helps deconfinement of bosonic spinons near the quantum critical point of 
the SO(5) WZW theory. We solve this effective field theory in the Eliashberg framework, and find 
non- Fermi liquid physics in thermodynamics and transport, where z = 3 criticality with dynamical 
exponent z plays an important role for hole dynamics. We discuss validity of our field theory, ap- 
plying it to a doped spin chain and comparing it with the slave-fermion framework. Furthermore, 
we discuss instability of the anomalous metallic phase against superconductivity and density waves 
of doped holes, resulting from competition between gauge and valance bond fluctuations. 
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I. INTRODUCTION 

Doping to an antifcrromagnetic Mott insulator has 
been one of the central interests in modern condensed 
matter physics, associated with high T c superconduc- 
tivity. Since the normal state is extremely anoma- 
lous, particulary shown from the absence of quasipar- 
ticle excitations near (±7r, 0) and (0, ±tt) momentum 
points [l[ and temperature quasi-linear behavior in trans- 
port experiments [2| although well defined electron ex- 
citations seem to exist near (±^,±§)[J, emergence 
of rather conventional BCS-type superconductivity from 
such an abnormal normal state has been a long stand- 
ing puzzle and still does. Such an anomalous behavior 
would be surely due to strong correlations between elec- 
trons, associated with Mott physics. 

Slave-boson approach has been one of the canonical 
frameworks for study of strongly correlated electrons. In 
particular, a doped Mott insulator problem was formu- 
lated in the slave-boson context, [3j where strong repul- 
sive interactions cause so called the single occupancy con- 
straint naturally imposed in the slave-boson representa- 
tion, and link variables arise as collective "order param- 
eter" excitations formulated as gauge fields. U(l) slave- 
boson gauge theory has been enjoyed both intensively 
and extensively for the doped Mott insulator problem, 
but such a formulation turns out to have fundamental 
difficulty for d-wave superconductivity emerging from a 
doped Mott insulator. [J] 

Wen and Lee have developed an SU(2) formulation, 
which extends the U(l) slave-boson theory to include 
fluctuations between nearly degenerate U(l) mean-field 
states, well applicable in underdoped regions. [1] Based on 
their SU(2) construction, they could obtain the doping 
independent decreasing ratio of superfluid weight with 



a confinement ansatz. In addition, they predicted a spe- 
cial structure of a vortex, which has a staggered flux core. 
Furthermore, they argued similarity between their stag- 
gered flux phase of the SU(2) slave-boson theory and the 
Gutzwillcr projected BCS wave function based on an ex- 
plicit numerical evaluation. 

Although the SU(2) slave-boson framework has ex- 
plained many kinds of aspects for high T c cuprates 
such as phase diagram, thermodynamics, transport, and 
etc.,[3j antiferromagnetic spin fluctuations are difficult to 
take into account in this context. High T c superconduc- 
tivity is sometimes argued to emerge from the spin liquid 
phase described by the spin sector of the slave-boson the- 
ory in the slave-boson community, but its connection to 
antifcrromagnetism is an important issue since the orig- 
inal problem is doping to an antiferromagnetic Mott in- 
sulator instead of doping to the spin liquid one. This 
motivates us to propose how to introduce antiferromag- 
netic correlations in the SU(2) slave-boson framework. 

Study of quantum antiferromagnetism associated with 
high T c cuprates has been performed in the context of 
spin liquid. In the bosonic representation of spin the 
half filled quantum antiferromagnet on the square lattice 
is described by the 0(3) nonlinear a model, allowing a 
quantum phase transition from an antiferromagnet to a 
quantum disordered paramagnet. Bernevig et al. have 
claimed that although the appropriate off-critical elemen- 
tary degrees of freedom are given by either spin 1 exci- 
tons (gapped paramagnons) in the quantum disordered 
paramagnet and spin 1 antiferromagnons in the antifer- 
romagnet, at the quantum critical point such excitations 
should break up into more elementary spin 1/2 excita- 
tions usually called spinons. 5] This was challenged by 
Senthil et al.j6| They argued that since the phase transi- 
tion in Ref. [5[ is supposed to fall into Landau-Ginzburg- 
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Wilson paradigm, the spectrum at the quantum critical 
point should be fully understandable only in terms of 
spin 1 bosonic degrees of freedom. 0] Senthil et al. pro- 
posed, as a possible candidate for a deconfined quan- 
tum critical point, a direct quantum phase transition be- 
tween a Neel antiferromagnet and a valance bond solid 
state, where one gets spinon condensation in the Neel 
state while instanton excitations (tunnelling events be- 
tween energetically degenerate but topologically inequiv- 
alent gauge vacua in the CP 1 representation of the 0(3) 
nonlinear a model) should possibly arise in the param- 
agnetic phase, whose condensation does not allow spinon 
dcconfincmcnt. Senthil et al. demonstrated that such 
proliferation of instantons is not the case at the quan- 
tum critical point, where a topological 9 term usually 
referred as a Berry phase term makes instantons irrele- 
vant and accordingly, makes it possible to achieve spinon 
dcconfincmcnt . 

This proposal motivated direct numerical simulations 
of various microscopic models to find such an exotic 
quantum critical point beyond the Landau-Ginzburg- 
Wilson paradigm. Actually, Sandvik has claimed that 
such a critical point exists indeed in a modified Hciscn- 
berg model with four-spin interactions. Q Furthermore, 
he pointed out an important thing, that is, only one 
length scale seems to exist and accordingly, critical ex- 
ponents for both staggered spin correlations and valance 
bond fluctuations are the same as each other, appar- 
ently in contrast with the original proposal for deconfined 
quantum criticality of the bosonic field theory. 

Resorting to the fermion representation of spin, QED3 
or QCD3 was obtained as an effective field theory at half 
filling, where its conformal invariant fixed point is al- 
lowed in the large N limit, identified with algebraic spin 
liquid. [§, [l(| An important point is that this effective 
field theory has an enlarged symmetry compared with 
its microscopic Hamiltonian. QCD3 resulting from the 
7r flux ansatz has Sp(4)«SO(5) while QED3 arising from 
the staggered flux gauge exhibits SU(4), where both are 
basically associated with spin and nodal structures. [9l[Toj 
This enlarged symmetry has important physical implica- 
tion that symmetry equivalent operators have their same 
scaling dimension, which should be taken into account for 
an effective field theory of such composite field variables. 

Actually, Tanaka and Hu have derived an effective field 
theory for competition between antiferromagnetism and 
valance bond solid. fll| Starting from the ir flux phase, 
they obtain QCD3 as its effective field theory although 
they do not take such gauge fluctuations into account 
explicitly. Their crucial observation is that the antiferro- 
magnetic order parameter is symmetry equivalent to the 
valance bond one via chiral rotation in SO (5). Introduc- 
ing a chiral rotated mass term with an SO (5) superspin 
vector, they could derive an SO (5) Wess-Zumino-Witten 
(WZW) theory. Validity of this description is further 
supported by the fact that the SO (5) WZW theory can be 
regarded as a natural extension of the SO (4) WZW the- 
ory for one dimension, [l2| which nonabelian bosonization 



for spin degrees of freedom or equivalently SU(2) chiral 
anomaly gives rise to. This discovery is meaningful since 
it tells a possible connection with the numerical simula- 
tion of Sandvik 8] although an appropriate treatment for 
the topological term is not known, thus true solutions 
for such an effective field theory are not found. Further- 
more, such competition between antiferromagnetism and 
valance bond solid is argued to be reflected as the checker 
board pattern in the scanning tunnelling microscopy. [l3j 

In this paper we follow the same strategy with the 
study of Tanaka and Hu for spin degrees of freedom, 
but away from half filling introducing charge fluctua- 
tions based on the SU(2) slave-boson framework. Two 
important problems can arise away from half filling. One 
cautious theorist may concern that a chemical potential 
term will appear in the Dirac theory of the fermionic spin 
representation away from half filling, thus breaking the 
original structure for the SO (5) WZW theory. However, 
such a term does not arise at least for the staggered flux 
phase of the SU(2) slave-boson theory since spinons are 
still at half filling even away from half filling. This is the 
reason why the SU(2) slave-boson framework is utilized 
in the present study. The other is how to construct cou- 
plings between doped holes and spin fluctuations. The 
second problem is our key issue. 

We find that it is not easy to derive such couplings if we 
resort to the nonlinear a model descri ptio n for the holon 
sector, derived in the previous study, [lj] In the present 
paper we perform fermionization for the holon sector via 
flux attachment, where a Dirac theory with a chemical 
potential term due to the presence of finite density of 
holons is found in the staggered flux ansatz. Then, we can 
construct couplings between spin fluctuations and holons 
in the same way as the spinon sector since its theoreti- 
cal structure is basically the same as that of the spinon 
part. Performing the standard gradient expansion, we 
have an effective field theory for SO (5) spin fluctuations, 
fermionic holons and their couplings. As a result, dynam- 
ics of doped holes in the antiferromagnetically correlated 
spin background turns out to be described by the SO (5) 
WZW theory for spin fluctuations, non-relativistic QED3 
around Dirac nodes for doped holes, and interactions be- 
tween valance bond fluctuations and holes. 

Fermionic holons remind us of the slave-fermion repre- 
sentation. We compare our effective field theory with the 
U(l) slave-fermion theory of the t-J model, and discuss 
similarities and differences between them. Then, we solve 
the effective field theory in the Eliashberg framework, 
where momentum dependence in the fermion self-energy 
and vertex corrections are neglected, justified by Migdal 
theorem and large N approximation with the fermion fla- 
vor number N at least for the second assumption. [l5[ The 
Eliashberg framework allows us to construct Luttinger- 
Ward functional of its closed form, thus making it pos- 
sible to solve the field theory self-consistently. Based on 
this framework, we discuss its phase diagram, thermo- 
dynamics, and transport near and away from the quan- 
tum critical point of our effective field theory, where non- 
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Fermi liquid physics is found, analogous with that of the 
slave- fermion framework proposed recently. [iH \vi\ 

An important feature of our field theory is the coupling 
term between valance bond fluctuations and doped holes. 
Considering that valance bond fluctuations are deeply re- 
lated with monopole excitations of staggered U(l) gauge 
fields in the bosonic field theory for spin fluctuations . [ 1 81] 
the presence of such a coupling term allows us to investi- 
gate how dynamics of doped holes affects deconfinement 
of bosonic spinons. Remarkably, we find that such cou- 
plings help spin fractionalization near the quantum crit- 
ical point of the SO (5) WZW theory. 

We also apply our field theory to one dimensional an- 
tiferromagnetic doped Mott insulator, where spin fluc- 
tuations are described by the SO (4) WZW theory while 
charge excitations are represented by relativistic QED2. 
Since one dimensional study is well known, possible con- 
clusions of our field theory in one dimension can justify 
validity of our description. Furthermore, reliable ana- 
lytic techniques can be utilized in one dimension, thus 
this study would suggest an important feature how doped 
holes affect spinon deconfinement. 

Finally, we discuss instability of the non-Fermi liquid 
metal phase against superconductivity and density waves 
of doped holes. As will be seen explicitly in the effective 
field theory [Eq. (10)], interactions between doped holes 
turn out to be mediated by both gauge and valance bond 
fluctuations. As a result, emergence of superconductivity 
from the non-Fermi liquid state is determined by compe- 
tition between gauge and valance bond interactions. For 
particle- hole channel instabilities, we discuss that the ho- 
mogeneous metallic phase can be stabilized at least when 
the flavor number of fermions (doped holes) is sufficiently 
large. 



II. EFFECTIVE FIELD THEORY 

Dynamics of doped holes in the antiferromagnetically 
correlated spin background is described by the t-J Hamil- 
tonian 

H = ~t I>I c j* + H.c.) + JJ2(S> ■ Sj - \n inj ). 



rewrite the t-J model in terms of these fractionalized ex- 
citations with hopping and pairing fluctuations 



(ij) 



(1) 



Introducing an SU(2) slave-boson representation for an 
electron field 

c iT = -L/4^. + = -j=(blja + by} 2 ), 

Cli = ^=h\iH- = - (2) 

where tpi+ = ^ J an d ^i- = ( ^ are SU(2) 
spinon-spinors and hi = ( J is holon-spinor, one can 



L = L + L S + L h , L = JrY, tT [ U ij U * 



13 h 



(ij) 



(3) 



i 

(ij) 

Lh = yih\(d T - fi- ia%T k )hi 

i 

+t r ^{h\u, j h J +H.c), 

(ij) 



where the SU(2) matrix field is C/j. 



and J r = and t r = | are redefined couplings. [|[ Since 
this decomposition representation enlarges the original 
electron Hilbert space, constraints are introduced via La- 
grange multiplier fields a* Q with k = 1,2,3. 

In the SU(2) formulation Wen and Lee choose the stag- 
gered flux gauge [3j 




U?f = -vV+?7 2 T3 exp[i(-l *r 3 ] (4) 

with a phase <f> = tan -1 ^^. Although the staggered 

flux ansatz breaks translational invariance, this formal 
symmetry breaking is restored via SU(2) fluctuations be- 
tween nearly degenerate U(I) mean-field states. For ex- 
ample, one possible U(I) ground state, the d-wave pair- 
ing one Uff c — ~x T 3 + (— f ) Iy+j!/ 77Ti can result from the 
SU(2) rotation Uff c = WiUf/w] with an SU(2) ma- 
trix Wi — exp|i(— f) lx+l! '| : ri|. Then, our starting point 
becomes the following effective Lagrangian 

i 

(ij) 

i 

■ i.y^iin ■ e. c .) 

(ij) 

+J r Ytr[U^Uf j p }, (5) 

where we have introduced only one kind of gauge field 
as important low energy fluctuations since other two 
ones, a* and are gapped due to Anderson-Higgs mech- 
anism in the staggered flux phase, [lij 

Our idea is to fermionize the holon sector attaching a 
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fictitious flux to a holon field 



L sf = 7] ^ dr ~ V ~ ia iO T 3)Vi 

i 

(ij) 

-i^2cio(r)lT 3 r)i - — (d x Cy - d. 



(6) 



where a bosonic field variable hi now becomes a fermionic 



one Tji 



Vi2 



with O 



It is important to notice 



that our flux attachment is performed in an opposite way 
for each isospin sector, confirmed by the presence of T3 
in 20(77] r 3Vi) = d x c y — d y c x . As a result, there is no net 
flux in the mean-field approximation of this construction, 
considering that the density of bn bosons is the same as 
that of bi2 bosons in the staggered flux phase. fl4j 

This observation is interesting since it suggests a con- 
nection with an SU(2) slave-fermion representation. If 



!?„■ is shifted to a?„- 



the Chcrn-Simons flux is trans- 



ferred to spinons, turning their statistics into bosons. 
Then, we have a bosonic spinon description with a 
fermionic holon, nothing but the slave-fermion represen- 
tation. 

Performing the continuum approximation for the long 
wave length and low energy limits, we find an effective 
field theory in terms of only fermionic variables 14 



(7) 



Dirac structure [lOj results from the staggered flux ansatz, 
where both ip and 77 are 8 component spinors and 

CT 3 
-£T 3 

(Tl \ , / (72 



Dirac gamma matrices are 70 = 



7i 



, , and 72—! 



a 3 , is slave-boson 



U(l) gauge field, where a finite bare gauge charge e 
is introduced. c M is Chern-Simons gauge field with its 
statistical angle O = it. It is important to under- 
stand that spinons are still at half filling even away 
from half filling in the SU(2) formulation. The single- 
occupancy constraint in the SU(2) representation is given 
by Infix + /a/»2 + - b] 2 b i2 = 1. Thus, if the con- 

dition of (b^bn) — {b\ 2 b i2 ) — § with hole concentration 

8 is satisfied, we see (fjifn + j\ 2 h'i) — 1> i- e -i spinons 
are at half filling. As a result, a chemical potential term 
does not arise in the spinon sector. Actually, this was 
demonstrated for the staggered flux phase in the mean- 
field analysis of the SU(2) slave-boson theory. [Tij On the 
other hand, a chemical potential term appears in the 
holon sector to form a Fermi pocket around the Dirac 
node. 



To find an effective field theory for low energy spin fluc- 
tuations, it is necessary to consider physical symmetry of 
the spinon sector. We know that this symmetry is closely 
connected with both spin and Dirac spaces. Since spin 
SU(2) symmetry is hidden in the present 8 component 
representation, we consider the redundant representation 



of 16 components with a Dirac spinor ip 



h 

-ft 



following Ref. 10]. Noting that the 

group space is composed of G = G Di rac ® G gauge <g> G» p, 
we see 15 generators associated with SU(4) svmmetrvflO| 
given by I ® I ® a, 73 (g) I ® a, 75 ® I <8 a, 27375 ® I <8 I, 
73 <8 1~3 ® I, 75 ® t 3 <g> I, and 27375 T3 (8 a, where 

73 = ( I ) and 75 = ^°^ 2 ^ 3 = i ( - 1 ) satis fy" 
ing [7^,73(5)]+ = 0. This implies that symmetry equiva- 
lent operators via SU(4) have the same strength for insta- 
bility, i.e., the same critical exponent for each correlation 
function. Actually, this was intensively discussed in Rcfs. 

Recently, C. Xu and S. Sachdev have claimed existence 
of a novel spin liquid fixed point, where such an SU(4) 
symmetry is broken down to SO(5).[19] At this fixed 
point most relevant spin fluctuations are Neel vector and 
valance bond fluctuations, and they compete with each 
other since 10 generators of 7®J®ct, 73®/®^, 75®/®ct, 
and 17375 8> I <8> I remain and such spin fluctuations are 
symmetry equivalent operators via chiral rotation. 

In this paper we also focus on competition between 
Neel and valance-bond fluctuations. Following Tanaka 
and Hu,[ll[ we introduce an SO(5) superspin vector v = 
(vi,V2,V3,v<i,V5), where the former three components 
form Neel vectors and the latter two ones represent x and 
y valance bond fluctuations, and consider the fermion 
mass term — m^(v ■ T)^ with T = (a x ,<T y ,<T z , 273,275) 
for the spinon sector. On the other hand, the holon mass 
term becomes — m v fj(ij3 W4 + 175^5) since it does not have 
spin. As a result, we find the following Lagrangian for 
symmetry "breaking" 



(8) 



The next task is to perform integration of Dirac 
spinons and expand the resulting logarithmic action 
for the superspin vector. Based on the gradient ex- 
pansion method, Tanaka and Hu have derived an 
SO (5) nonlinear a model with a WZW term ignor- 
ing gauge fluctuations, [1JJ S sp m = SnlsM + Swzw, 
where Snl s m = J d 3 x-^(d fJ ,v k ) 2 and S WZ w = 

i Area(S±) So dt I ^ ' x ^abcdeV a dtVbd T V c d x VddyV e with 

Area(S 4 ) = £(5/2) • Since Dirac fermions are massive 
in the symmetry "broken" phase and their fluctuations 



5 



are ignored in the low energy limit, gauge-fluctuation 
corrections will be irrelevant considering that they can 
appear through fermion bubbles. Only one point that 
should be careful is a topological contribution, associated 
with an imaginary term. Although the WZW term is 
nicely derived in the absence of gauge fluctuations, an 
additional imaginary term may arise, a coupling term 
between a topologically nontrivial fermionic current 
and gauge field. [20j If we represent the Dirac spinor as 



4>n 



where \n is a 2 component spinor with 



an isospin index n = 1,2, one can see that each sector in 
the Dirac space gives rise to such a term. However, their 
signs are opposite, thus such terms are cancelled. [2l| 
This is well known to be cancellation of parity anomaly 
in condensed matter physics. Another way to say this is 
that the signs of mass terms for Dirac fermions (%+ and 
X^) are opposite, resulting in cancellation of the parity 
anomaly. 

Based on the above discussion, we reach an effective 
field theory 



S = J d 6 x 

d 3 x 



^9 



+ S- 



dimension. Second, dynamics of doped holes is described 
by U(l) gauge theory with finite density of fermionic 
holons around four Dirac nodes. Thus, non-Fermi liq- 
uid physics is expected naturally. Third, interactions be- 
tween spin fluctuations and doped holes emerge as cou- 
plings between valance bond fluctuations and fermionic 
holons. Considering that valance bond fluctuations are 
deeply connected with monopole excitations of CP 1 or 
staggered U(l) gauge fields in the bosonic field theory 
for spin fluctuations, this coupling form implies how dy- 
namics of doped holes affects spin fractionalization, i.e., 
deconfincment of bosonic spinons. This will be discussed 
more deeply. 



III. CONNECTION WITH U(l) 
SLAVE-FERMION THEORY 



Although the effective field theory Eq. (10) has the 
similar spirit with the slave-fermion framework, it has 
an important different point. To understand this more 
clearly, we consider the U(l) slave-fermion representation 



wzw 



m^iJ. - ia l T 3 - *c M r 3 - iA^)n - Hh'njo'n 

i 1 



(9) 



where spin fluctuations are described by the SO (5) WZW 
theory. An interesting observation is that the Chern- 
Simons contribution becomes irrelevant if the holon dy- 
namics is in a critical phase. Shifting the slave-boson 
gauge field as a 3 ^ — c M and performing integration of 
Chern-Simons gauge fields, we obtain ~ (d x d x a 3 ) ■ 
(d x a 3 ). This contribution is irrelevant since it has a 
high scaling dimension owing to the presence of an addi- 
tional derivative. Considering that the density of holons 
is finite to allow a Fermi surface (pocket around the Dirac 
point), it is natural to assume that the fermion sector is 
in criticality. We note that this kind of argument was 
well utilized previously. [22] As a result, we find an ef- 
fective field theory for an antiferromagnetic doped Mott 
insulator problem 



S= d 



^9 



+ / d 3 x 



+ 2 e 2 ( € ^i^ a j) 



wzw- 



(10) 



Several remarks are in order. First, the spin sector is 
described by the SO (5) WZW theory even away from half 
filling, starting from the SU(2) slave-boson theory in the 
staggered flux gauge. Validity of this description will be 
further supported, comparing the present effective theory 
with the slave-fermion framework and applying it to one 



(11) 



where ipi and bi a are fermionic holon and bosonic spinon, 
respectively. Inserting this representation into Eq. (1), 
we obtain the following expression for both exchange- 
interaction and electron-hopping terms, 

J ^(5, • §j - -mrij) 

ij 

- ^ |A>/ - J^^eaMf, + H.c), 

ij ij 
ij ij 



(12) 



where each composite field is given by 



a/3 a 

representing short-range antiferromagnetic correlations, 
ferromagnetic spin fluctuations, and hopping of holons, 
respectively. 

In the low energy and long wave-length limits one can 
set the above collective fields as 

A 4 , = Ae** , X % = Xbe ia » , x% = X^ ,(13) 

where amplitude fluctuations are frozen to be their 
saddle-point values, and only phase fluctuations are 
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kept importantly. Then, the resulting slave-fermion La- 
grangian becomes 

Leff = L s + Lh + L c + L , 

L S = J2 blidr - tik* - t Xb Yj>L'" : b jo + H.c.) 

i ij 

- JA^{e- lc ^e af3 b la b jfi + H.c), 

ij 

i ij 

L c = iY,^(blb ta +4^-2S), 

i 

L a = N L (JA 2 +4t X ^ X b), (14) 

where Xi in L c is a Lagrange multiplier field to impose 
the single occupancy constraint with S = 1/2, and Nl 
in Lq is number of lattice sites. 

Comparing Eq. (10) with Eq. (14), one will find sev- 
eral different points. First, the spectrum of holon exci- 
tations is not relativistic in Eq. (14). However, this can 
be adjusted, allowing ir flux in the U(l) gauge field ay, 
i.e., dij = 7T. Then, four Dirac nodes arise, and hole 
doping gives finite density of holons, resulting in the hole 
pockets around the nodes. 

Second, spin dynamics is described by usual spin 1 ex- 
citations in Eq. (10) while it is expressed with spin 1/2 
fractionalized spinons in Eq. (14). However, the pres- 
ence of the WZW term in Eq. (10) may give rise to 
spin fractionalization as one dimensional physics [12| or 
the previous proposal [6| for deconfined quantum critical- 
ity. If we start from Eq. (14), we focus on two kinds 
of gauge fluctuations, corresponding to staggered U(l) 
gauge fields cy and uniform U(l) gauge fields ay, re- 
spectively. In particular, the staggered U(l) gauge field 
cy turns out to be the same as the CP 1 gauge field if the 
Schwinger-boson effective Lagrangian is mapped onto the 
CP 1 gauge theory of the 0(3) nonlinear a model in the 
long wave-length limit. [18] In this respect, if we assume 
that staggered gauge fluctuations mediate confining in- 
teractions between spinons, we expect to see spin 1 fluc- 
tuations described by a er model-type theory. However, 
the fate of valance bond fluctuations is not clear in this 
case. 

Third, the coupling structure between spin fluctuations 
and doped holes in Eq. (10) differs from that in Eq. (14). 
In the slave-fermion Lagrangian Eq. (14) direct couplings 
between spin and charge degrees of freedom do not ex- 
ist although nonlocal current-current interactions appear 
in the low energy limit, associated with the single occu- 
pancy constraint. Even if we assume staggered gauge 
fluctuations confining, it is difficult to find a coupling 
term between spin fluctuations and doped holes. On the 
other hand, our effective field theory Eq. (10) exhibits its 
direct coupling term explicitly. This coupling term is an 
important ingredient for dynamics of doped holes in the 
quantum antiferromagnet. As demonstrated intensively 
in the previous work,[l8j valance bond fluctuations are 



deeply related with instanton (monopole) excitations of 
staggered U(l) gauge fields cy. Actually, va and v$ fields 
can be identified with an instanton operator in the pres- 
ence of a topological 9 term. This implies that the pres- 
ence of doped holes affects spinon confinement directly 
via their couplings. In the cumulant expansion for the 
coupling term one can see that holon dynamics modifies 
monopole-monopole correlations, giving rise to dissipa- 
tion. As a result, an effective theory for monopole exci- 
tations away from half filling is changed seriously from 
that at half filling. In the next section we will see that the 
presence of doped holes helps bosonic spinons deconfined 
near the quantum critical point. 

IV. PHYSICAL PROPERTIES 
A. Eliashberg framework 

Three kinds of field variables, that is, fermionic holon, 
SO(5) superspin vector, and uniform U(l) gauge field 
make our effective field theory complicated. In this re- 
spect it is not easy to treat such all degrees of freedom 
self-consistently. Recently, it was explicitly demonstrated 
that Eliashberg framework is the minimal self-consistent 
treatment for an effective field theory near its quantum 
critical point. [IB] The Eliashberg treatment neglects mo- 
mentum dependence of a fermion self-energy and vertex 
corrections. The first assumption is based on the fact 
that momentum dependence of a fermion self-energy is 
regular, and singular physics arises from its frequency de- 
pendence. This can be checked explicitly at least in the 
one loop level. The second assumption is more serious 
than the first one, sometimes called Migdal theorem. [23j 
When fermions are much faster than bosons, vertex cor- 
rections can be neglected since the pre-factor in the renor- 
malized vertex is given by the ratio of fermion and boson 
velocities. However, this turns out to be not sufficient 
for the Eliashberg framework. [l5| Another parameter is 
shown to be need, that is, the fermion flavor number N. 
In the large N limit the Eliashberg framework is justified. 

For the Eliashberg treatment we rewrite the effective 
field theory as follows 

S = S 7 ] -+- S v -(- S a -\- Si n t, 

S v = J d 3 x|^7 P (5p - iAf^r] - ^fnwj, 

S v = I d 3 x{j-(d^v k ) 2 +m 2 v (\v k \ 2 -l)}+S wzw , 

S a = j d 3 x[^ (e^a*) 2 }, 

Sint = /d 3 a;|-m^f?(i73«4 + i^v^r) - ia^-f^r]^, 

(15) 

where m 2 represents mass of superspin vector bosons 
arising from the SO(5) rotor constraint 53fc=i \ v k\ 2 = L 



7 



It will be determined self-consistently. 

Performing the cumulant expansion for Si nt , we find an 

effective action S eff = S v + S v + S a - |((S? nt ) - (S^t) 2 ). 
Although this expression is for the second order, it in- 
cludes an infinite order actually. Accordingly, we can 
construct the corresponding Luttinger-Ward functional 
in the Eliashberg framework 



F, 



lw — F LW + Fl w + Fl w + Y v 



d d k 
(2n) d 

-E 7) (iw)G r) (fc,icj) , 



tr 



In 



id J * ' m,n=l 
5 



in 



-U a (q,m)D a (q,in) 



Y, = 



) ^ r d d k ^ r d d q 



E tr[D™ n (q, in) lm G v (k + q,iui + iil)j n G v (k, iu)], 



a 2 if J it J ^ y 



(2n) d j^J (2tt)» 
tr[D% v (q, iil)^T 3 G ri {k + q.iuo + iO) 7l/ r 3 G^(fc, *w)]<16) 
Grj(k, iu>) in F^ w is a renormalized propagator for holons, 
given by G v {k,iLu) = ^g^ 1 (k, iu>) + S^(io;)| , where 

g v (k,iuj) = (i~/oU) + «7ifci + fJ-h^o^j is its bare prop- 
agator, and S^(iw) is its momentum-independent self- 
energy. D™ n (q,i£l) in F£ w is a renormalized propaga- 
tor for superspin vector fields, given by D™ n (q, iSl) = 

{d^fotf^mn + n^™(g,iO)} , where d„(g,ifi) = 

( ? 2g° + m ?>) i s its bare propagator, and n m ™(q, iQ) is 
its self-energy. D a (q, ift) in is a renormalized kernel 
for the gauge propagator D^ l '(q,iil) — D a (q, i£L)\Snu — 

Sf^j, given by D a (q,i(l) = {d^iq, iQ) + U a (q, 

where d a (q,iQ) = ( q ^ j is its bare kernel, 
and n a (<7, ifi) is its self-energy in IL£ v (q,iQ,) = 
H a (q, id) {Sfnv — ^r~^j- Y v is introduced for self-energy 
corrections resulting from the first term in Si nt of Eq. 
(15) while Y a is for those arising from the second term 
of Sint- In the self-energy functional Y v 74 should be 



replaced with 73, where this problem appears from our 
notation form. 

It is important to notice that Luttinger-Ward func- 
tional is not usually written in a closed form. However, 
the Luttinger-Ward functional can be written in its closed 
form at least for the Eliashberg framework. [24| Actu- 
ally, performing variation for the Luttinger-Ward func- 
tional with respect to each self-energy, i.e., S S ^ L ,) V .\ = 

Eliashberg equations 



5 Flw — — 0, we find self-consistent 



= m ? T E [ 4z 



d d q_ A 



d d q 

^ J (2n) d 



D™ n (q, in) lm G v {k F + q,iLU + itl) ln + T 
D^ v (q,iCl)j l± T3G ri (kF + q,iu) + ifl)j„T 3 , 

nr(^) = TW|^L ± 

iu m 7 n=4: 

( — ^tr [j m Gr, (k + q,iuj + iU)~f n G n (k, iu))]) , 

(-^tr [7^3(7,, (A: + q, iu) + ffl)7 !/ r 3 G n (fc, iu))f) . (17) 

The holon self-energy results from both valance bond and 
gauge fluctuations, where 74 — > 73 is performed. The su- 
perspin vector self-energy arises from holon fluctuations, 
where 74 is also replaced with 73. Notice II™™ (q, iQ) = 
for m, n = 1,2,3. The gauge field self-energy appears 
from holon current fluctuations. Eqs. (16) and (17) com- 
plete the Eliashberg framework. 

B. Simplification of Luttinger-Ward functional 

Using the Eliashberg equations (17), one can simplify 
the Luttinger-Ward functional Eq. (16) as follows 



LW 



= - T E 



d d k 

(2n) d 



tr 



In 



Y, T) (iu))G ri (k, iu) 

T E / t|^[ E ln{Cfoi«)*m* + lC B (<M«)} 
iQ ^ \ ' m,n=l 

W ' + T E/ (0 ln { rf a 1 (9,*^) + n a ( (Z , l ^)}, (18) 



where the self-energy parts for superspin vector and U(l) 
gauge fields are cancelled out from the use of Eq. (17). 
Then, one can see that the holon free energy is nothing 
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but the free energy of Fermi liquid as follows 

d d k 



N vPri Tj2\u\-- (j 



(19) 



the final expression for the Eliashberg free energy 

^N, lPvrr2 C~ 2 T >p f d d q 
2 9 + |f 7 (2tt)<* 

|fi| 
'7 



{3 ln(g 2 + ft 2 + r 2 ) + 2 ln(V + £~ 2 + 7„ 



Here, /?,, is the density of states around the Dirac node, 
and the number of Dirac nodes is N v = 4. f(y) 



T ^ [ d i , ( 2^ M 



)} 

(24) 



IS 



the Fermi-Dirac distribution function with temperature 
scaling. As a result, we find the Eliashberg free energy 
of our effective field theory 



where 



r 



ngN^m^Pr, 'ne 2 N rj p ri 
2gm v , j v = , 7 a = 



Flw = 



* N VPV T 2 2 

6 v 

jd 



+ T E / E ^{d-\ q ,in)s mn + ur(q,in)} 

in J \ ' m,n=l 



(20) 



The remaining thing is to evaluate each self- 
energy. The superspin vector self-energy H™ n (q,iQ) = 
Tl v (q, ity&mn for to, n — 4, 5 is found to be 



U v (q, in) 



4 ^ 



(21) 



This is nothing but the standard Landau damping term, 
originating from fermion excitations near the Fermi sur- 
face. Since the density of holons is finite due to hole 
doping, emergence of the Landau damping term is quite 
natural. Then, the renormalized propagator for super- 
spin fluctuations is given by 



for m, n = 1, 2, 3, 



represent the correlation length for superspin fluctua- 
tions, Landau damping coefficient for superspin fields, 
and that for gauge fields. 

Several remarks about the Eliashberg framework are 
in order. First, the Landau damping contribution does 
not depend on the fermion self-energy in the Eliashberg 
framework, allowing one to use a bare fermion propaga- 
tor instead of its full green's function. [l5| Second, z — 3 
criticality supports Migdal theorem since on-shell (or res- 
onance) fermion momenta are larger than on-shell boson 
one at the same energy. Third, the number of Dirac nodes 
Nfj plays the same role as the fermion flavor number in 
the Eliashberg framework. In this respect the Eliashberg 
framework works well for our effective field theory. 



C. Phase diagram 

Performing variation of the free energy Eq. (24) with 
respect to the correlation length, i.e., = 0, we ob- 

tain the self-consistent equation for the correlation length 
in the Eliashberg framework 



2y 



2 



n 2 



D™{q,itt) 



q 2 + £-2 + lv 



(25) 



2ff 



nN^m^fri \n\ 



for to, n = 4, 5, (22) 



implying that antiferromagnetic spin fluctuations are de- 
scribed by z = 1 theory while valance bond fluctuations 
are expressed by z — 3 theory, where z is the dynami- 
cal exponent. The gauge self-energy is also given by the 
Landau damping term H a (q : iQ) = vN ^ Pt ' M ; thus the 
renormalized gauge propagator becomes 

1 



Notice that interactions between valance bond fluctua- 
tions and holons result in the z — 3 part. 

Performing integration, we find the following expres- 
sion for the correlation length 



A \ 



^Tln{sinh( ^ r ^ } 
3tt I V 2T // 



3tt 
-2T 



^- a -^T-(^ln(l- e -^ 



4T 



D a {q,iQ) 



q 2 TiNnPr, \S1\ > 
2(? 4 q 



(23) 



-(^7„)- 1 ln{sinh( 



(C 3 7.)" 



2T 



}+^E 



OO _ l. (| 

e t 



fc=i 



k 2 



described by z = 3 critical theory. 

Inserting bosonic self-energies into Eq. (20), we find 



2.9 



T 



In 



sinh^ 



2t)- 1iIMh1j 



(26) 
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AF&VBS y''' T ~^/2 
z=l&z=3 / VBS 
/' B z=3 



CO QCP QD 



A 



FIG. 1: (Color online) Schematic phase diagram : CO, QD, 
and QCP represent classically ordered, quantum disordered, 
and quantum critical point, respectively. Two crossover scales 
emerge as T 



and T ~ 2~ j identifying three 

—n ^ both z = 1 

AF and 2 = 3 VBS fluctuations are gapped while in regime 



regimes. In regime (A) of T < ^ "'^ — < - 



(C) of 



< T both spin fluctuations are criti- 



cal, that is, in the quantum critical regime. In regime (B) of 

^ — < T < only 2 = 3 VBS fluctuations are critical, 
and 2 = 1 AF ones are gapped. 



As shown in this equation, there exist three regimes, (A) 



T < 



< 



; (B) 



< T < 



-, and (C) 



< < T, emerging from coexistence of z = 1 
(antiferromagnetic) and z = 3 (valance bond) fluctua- 
tions. In regime (A) both 2 = 1 and z = 3 fluctuations 
are gapped while in regime (C) both spin fluctuations are 
critical, that is, in the quantum critical regime. In regime 
(B) only valance bond fluctuations (z = 3) are critical, 
and z = 1 antiferromagnetic ones are gapped. The phase 
diagram is shown in Fig. 1. 



Considering 



near the quantum critical 



2^ ^2 

point, we can simplify the above equation, and find that 
antiferromagnetic fluctuations play an important role in 
determining the SO (5) superspin correlation length near 
the quantum critical point. As a result, we obtain the 
following expression for the correlation length 



r\ 9 



9c;T<'—) 



{( 



i 



2tt\ 



9 J 

3 \g c g 



A 



2tt -i 
Yg) 



r 1 



(27) 



where g c = \^ + ^ )A is the quantum critical point. [251] 

Recalling g(S) oc m 2 (5) in the gradient expansion and 
m 2 (5) (x \S — 5q\ with hole concentration S < 5q near the 
quantum critical point, we find the critical hole concen- 
tration in the Eliashberg framework 



So 



|{(S + 5W- 



(28) 



where c is a positive numerical constant. Thus, an anti- 
ferromagnetically ordered phase appears in 5 < S c , and 
a quantum disordered state arises in 5 > 6 C . The na- 
ture of the disordered state will be determined by SO (5) 
symmetry breaking effective potentials. 



D. Deconfined quantum critical point 

We consider how interactions between valance bond 
fluctuations and holons affect the deconfined quantum 
critical point, proposed to exist at half filling. As dis- 
cussed before, valance bond fluctuations can be identi- 
fied with monopole excitations of staggered gauge fields. 
Introducing ^ = -^(i^+zt^) and considering SO(5) sym- 
metry breaking for the WZ W term as follows in Ref . [ll[ , 
we find an effective field theory for monopole or valance 
bond fluctuations from the effective field theory Eq. (10), 



irN^m^pr, |fi| 



2fl 



*(?,»n) 



J dr J d 2 ry m (^ i + ¥ i ). 



(29) 



Here, the cubic power in the last term results from the 
WZW term with SO (5) symmetry breaking, where y m is 
the monopole fugacity. If the topological 9 term is not 
taken into account, the condensation-induced term will 
be given by - $ dr J d 2 ry m (* + ft). 

An important point is that dynamics of valance bond 
excitations is described by z = 3 critical theory at the 
quantum critical point. As a result, two spacial di- 
mensions are already in the upper critical dimension, 
thus higher order interactions beyond the Gaussian term 
are irrelevant. It means that the WZW-induced cubic 
term can be neglected safely at the quantum critical 
point. Equivalcntly, the monopole fugacity vanishes at 
the quantum critical point, indicating deconfinement of 
bosonic spinons. We point out that the topological term 
plays an important role for deconfinement even away 
from half filling. If the topological term is ignored, the 
monopole- fugacity term of linear in ^ is relevant, giving 
rise to confinement. 

Although proliferation of monopole excitations is pro- 
hibited due to fermion excitations, monopoles arc in a 
critical state of z = 3 instead of a gapped phase. Such 
critical valance bond fluctuations will affect holon dy- 
namics. This feedback effect to holon dynamics is an 
important ingredient for the source of non-Fermi liquid 
physics near the quantum critical point. 
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E. Thermodynamics 



F. Transport 



We study thermodynamics in each regime. We start 
from the following expression for the Eliashberg free en- 
ergy 



We discuss electrical transport of the effective field the- 
ory Eq. (10). As shown in the Eliashberg equations (17), 
the holon self-energy is given by two contributions 



ttN vPv 2 £ 2 
F LW = —T — 

*"■{"■* (^F)} 

1 Z" 00 y \ POO 

-2^1 d " COth (^)j Q *?t» _1 (7. ? ). (30) 



The first term is due to the Fermi liquid contribution 
of holon excitations, and the second is the correlation 
length term. The third contribution results from z — 
1 antiferromagnetic fluctuations while the fourth term 
arises from z = 3 valance bond excitations. The last 
term originates from z — 3 gauge fluctuations. 

Considering the specific heat coefficient given by 

j(T) = ^P- — — 9 F g^ 7 " > , we obtain its analytic ex- 
pression in each regime, 



l{T > n ) 



>8/3 



V2vr 



+ s)(^ /3 



2/3 
la' 



dy{- 



u <- sinh y 



y 5/3 + y 8/3 cothy 



7 (T< 



sinh y 
1 



T -l/3 



V27T 2 + 8W 7a 



y5/3 y8/3 CQ^y 



sinh y sinh y 



r -l/3 



(31) 



where only dominant contributions are shown. As dis- 
cussed previously, we have three regimes, (A) T < 



where both antiferromagnetic and valance bond 



< T < where 



fluctuations are gapped, (B) ^ 1 ^ — 

only valance bond fluctuations are critical, and (C) < 
T where both antiferromagnetic and valance bond fluctu- 
ations are critical. In the regime (A) contributions from 
superspin fluctuations exhibit an exponential dependence 
of temperature and ignored in the low energy limit. Dom- 
inant contributions are driven by z = 3 critical gauge 
fluctuations, resulting in ~f(T) ~ T -1 / 3 . In the regime 
(B) antiferromagnetic fluctuations cause an exponential 
dependence of temperature while both valance bond and 
gauge fluctuations give rise to j(T) ~ T" 1 / 3 due to 
their z = 3 criticality. In the regime (C) z = 3 criti- 
cal valance bond excitations and gauge fluctuations allow 
7(T) ~ T" 1 / 3 while z = 1 critical antiferromagnetic fluc- 
tuations result in j A f{T > = § f/ °° dx-^r^T, 
sub- leading thus ignored in the low energy limit. 



5 

E 

m.n—A 



D™ n (q, ift)-/ m G v (k + q,ioj + in)j n , 



m J v ; 

D^(q, ifl) ltl T 3 G v (k + q,iuj + itt)y u T 3 , (32) 

where the first is due to valance bond or monopole fluc- 
tuations, and the second comes from gauge fluctuations. 

We first consider the self-energy due to valance bond 
fluctuations. The above expression can be written as 
follows 



Mfc <£J (27T) 



d d q 



q2 _|_ £-2 _|_ iu) + iQ — qcost 



(33) 



where the self-energy correction of doped holes via z = 3 
valance bond fluctuations is clearly seen. Here, kp = 



k^ 2 + ky 2 — fih is the holon Fermi momentum. Then, 
the imaginary part of the self-energy is given by 



gm 2 

27T 3 /i, l 



+ ijikf + /z/,70] 



M 



<1 



C?fil / dq : nr. 



where Wick rotation is performed at zero temperature in 
order to see frequency dependence of the self-energy. 

Performing momentum and frequency integrals, we 
find 



2 ' 4V3v 7 y 3 



9 m ri I |2/3 

7oM 1 > 



3£> < 



-) 



Note that the |w| 2//3 behavior is the hallmark of z = 3 
criticality in two dimensions. [2^] The self-energy correc- 
tion due to gauge fluctuations also gives rise to Sy£"(u;) cx 

M 2 / 3 . 

At finite temperatures the zero-frequency self-energy 
corrections turn out to diverge in the one-loop ap- 
proximation. However, such divergences due to both 
gauge and valance bond fluctuations need not be given 
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much attention because such self-energies are not gauge- 
invariant, thus they do not have any physical meaning. 
These divergences should be considered as an artifact of 
gauge non-invariance. Gauge invariance can be incorpo- 
rated via vertex corrections, which cancel the divergent 
parts in the self-energie s, g iving rise to gauge invariant 
finite contributions. [27l l28t l29j This corresponds to the 
transport time, given by q 2 ~ T~ multiplication in the 
quasiparticle life time. As a result, we find the following 
expression for the electrical resistivity 



p(T) <x T 4 / 3 , 



(35) 



consistent with the previous results. [261] We obtain non- 
Fermi liquid physics in both quantum critical and disor- 
dered phases, where both valance bond and gauge fluctu- 
ations cause the non-Fermi liquid transport in the quan- 
tum critical regime while only gauge fluctuations result 
in that in the disordered phase. 



APPLICATION TO ONE DIMENSION 



Although our effective field theory Eq. (10) is physi- 
cally reasonable, we would like to justify its validity ap- 
plying it to one dimension. In one dimension the spin 
sector is described by the SO (4) WZW theory, and the 
charge sector is represented by QED2 without the chem- 
ical potential term. Accordingly, the coupling term be- 
tween valance bond fluctuations and holons is adjusted. 
The resulting effective field theory in one dimension is 
obtained to be 



2tt 



Area(S 3 ) J 
d 2 x 



fc=i 
1 

dte abcd v a dtv b d T v c d x Vd 



1 



(36) 



where 7 M and 75 are 2x2 matrices, and rj is a four com- 
ponent Dirac spinor. The SO (4) WZW theory of the spin 
sector has been derived in Ref. [II], using the path in- 
tegral formulation for non-abelian chiral anomaly. Here, 
we investigate the role of massless Dirac fermions, intro- 
duced from Eq. (10) by reducing the two dimensional 
theory down to one dimension. 

Performing the abelian bosonization for the fermion 



sector, we obtain the following expression 



s=j d 2 4^x>^) 



2tt 



Area(S 3 ) J 
d 2 



k=l 
1 

dte abcd v a dtV b d T v c d x Vd 



+ (— ViSm^V^Trifi+J + (^— m r) ^ va sin^V^^-^ 



1 



■iA^-^e^ud^ 



1 
1 



3\2 



(37) 



where the subscript ± in the bosonic field <p± represents 
the SU(2) doublet of T3, and A is a cutoff associated with 
band linearization. 

Performing integration for U(l) gauge fields, we find 
a mass-type term -^{<f) + — </>-) 2 . This allows us to set 
(/>+ = (j>— = <j) in the low energy limit. Shifting y/Ancj) 
with — f + y/^rO, we are led to 



r 1 4 
S = J d 2 x[-J2(d,v k ) 



+i 



2tt 



Area{S 3 ) J 



1 

dte abcd v a dtv b d T v c d x Vd 
'2A 



-iA^y—efiudvO 



(38) 



It is interesting to see that valance bond excitations lead 
to charge density wave fluctuations, consistent with our 
expectation. 

The valance bond and charge density wave coupling 
term can be taken into account in the cumulant expan- 
sion, and the correction part is given by 



8S=- 



1/2A \2 

= -2lT m " 



fx / d 2 x' 



{t> 4 (x)<cosfV47r0(a:)) cos^Vl^a/)))^') 

+ cos(yiTT9(x)^ (v4,(x)va(x')) cos (V4^(.t'))} 
= SS Vi + 6S , (39) 



where Si nt = — J d 2 xy^-m ri jv4 : cosy\/ATr9j is the cou- 
pling term. 

It is not difficult to evaluate the density-density corre- 
lation function since charge fluctuations are described by 
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the noninteracting Gaussian ensemble if metallic charge 
dynamics is assumed. In this case we find 



spin liquid Mott insulator is expected to appear in this 
case. 



4tt6(x)) cos(V4k6(x'))) cx cosh(47r(0(x)0(a;')> 



cosh(47rC e ln|:r-.cc'|) -> \x - x'\ inCe , 



(40) 



where Cg is a positive numerical constant, and the last 
part is valid at large distances, i.e., \x — x'\ — ► oo. 

Inserting this expression into the spin sector, we obtain 
an effective theory for SO (4) spin fluctuations 



Si,. — I d x 



+ i 



y fc=i 



2?r 



Area(S 3 ) J 



dt€abcdVad t V b d T V c d x V d 

-2A 



(41) 



where C„ 4 is a positive numerical constant. An impor- 
tant point is that metallic charge fluctuations give rise 
to confining interactions between monopoles, suppressing 
monopole fluctuations. This tendency is completely con- 
sistent with our previous two dimensional analysis, where 
holon fluctuations cause dissipative monopole dynamics 
described by z — 3, prohibiting their proliferation. In one 
dimension charge dynamics suppresses monopole conden- 
sation more strongly. 

An immediate question is the nature of spin dynamics 
descried by Eq. (41). Since monopole excitations will 
be suppressed via charge fluctuations, spinon deconfinc- 
ment is expected to appear. A relevant point is whether 
spin fluctuations are critical or not. The SO(4) WZW 
theory is well known to exhibit criticality without the 
fcrmion-induced monopole-suppressing term at half fill- 
ing, where monopole excitations turn out to be irrelevant 
due to the presence of the topological WZW term, which 
differs completely from the charge-dynamics monopole- 
suppressing mechanism. Since the SO (4) symmetry is 
broken by the presence of charge fluctuations shown by 
the last term of Eq. (41), we expect that the WZW 
term may not be relevant in this case, thus resulting in 
spin gap. This seems to be consistent with our physical 
intuition that charge fluctuations will cut spin correla- 
tions, making their correlation length short. Then, the 
resulting state is identified with the Luther-Emery phase, 
where spin fluctuations are gapped and charge excita- 
tions exhibit enhanced superconducting correlations. [30j 

On the other hand, if charge fluctuations are gapped, 
i.e., in the Mott insulating phase, their density-density 
correlations will vanish at large distances as follows, 

(cos where 

£~ is associated with their excitation gap. Then, spin 
dynamics will be described by the pure SO (4) WZW the- 
ory in the long wave-length limit. As a result, the critical 



VI. SUPERCONDUCTIVITY AND STABILITY 
OF THE ANOMALOUS METAL 



Until now, we have seen the nature of the anomalous 
metallic phase. An important problem is how d-wave su- 
perconductivity evolves from the non-Fermi liquid state. 
In our effective theory approach superconductivity can 
appear only from pairing of fermionic doped holes. As 
shown explicitly in our effective action Eq. (10), interac- 
tions between doped holes are mediated by both gauge 
and valance bond fluctuations. An important thing is 
that these interactions compete with each other. 

Since holons carry an isospin quantum number repre- 
sented by the T3 matrix, two kinds of fermion pairs can be 
considered. When their isospin quantum numbers are dif- 
ferent from each other, gauge fluctuations cause attrac- 
tive interactions for the particle-particle channel while 
valance bond fluctuations give rise to repulsive ones. This 
is nothing but the pairing possibility between b\ and 62 
bosons in the bosonic description [14j of the SU(2) slave- 
boson theory, discussed by Wen and Lee.Q In this re- 
spect such pairing possibility is the unique feature of the 
SU(2) slave-boson description for superconductivity. 

On the other hand, when their isospin quantum num- 
bers are the same as each other, gauge fluctuations 
give rise to repulsive interactions for the particle-particle 
channel while valance bond fluctuations cause attractive 
ones between holons on different sublattices. Do not con- 
fuse the present uniform (ferromagnetic) gauge fluctua- 
tions with staggered (antiferromagnetic) ones. Staggered 
gauge fluctuations, discussed previously in the slave- 
fermion context, cause attractive interactions between 
doped holes on different sublattices. Gauge fluctuations 
appearing in the present theory are in the uniform chan- 
nel associated with the third component of the SU(2) 
slave-boson theory. The presence of T3 matrix in the 
gauge coupling supports this fact. As a result, such 
holons have the same gauge charges, repulsing each other. 

Based on the above discussion, we conclude that emer- 
gence of superconductivity is determined by the rela- 
tive strength between the couplings of gauge and valance 
bond fluctuations with doped holes. An interesting ob- 
servation is that the presence of repulsive interactions in 
both cases may favor d-wave pairing of doped holes. 

Now, we derive a full set of Eliashberg equations in- 
cluding the pairing channel. In this paper we consider 
only the first case, that is, pairing between different 
isospins. One can construct the Luttinger-Ward func- 
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tional in the same way as the previous case Eq. (16), 
Flw = F2 W + Fl w + F£ w + Y v + Y a , 



d d k 

-S,(iw)G^(fc,iw)l , 
/ d<iq 



ti- 



ll! 



in J ' m,n=l 

5 



m.n— 1 



m 2 ,, 



-n o (g,iO)£) a ( g ,ifi)], 

2 ^-7 (^Y (^) d 



£ tr[LC n ( g , iO) 7m G 7) (fe + g, l ^ + infrnG^A;, iw)], 

m,n— 4 

a 2 trV w ^7 w 

tr[D^(q, ifl)j^T 3 G v (k + q, iuj + ifi) 7 (fc,iw)](42) 

Here, G,j = g _1 (fc, 2<j) + E(fe, iw) is the ma- 

trix green's function with an anomalous propagator, 
where g(fc, iu/) = g(k,iu)l with g (k,iu>) = {i^f^k^ + 
MftTo) -1 is the bare green's function and E(fc, iu/) = 
T,(k,iuj) A(k,iu) 
A*(k,iuj) Ti(k,iu>) 
anomalous part. 

The above Luttinger-Ward functional results in the fol- 
lowing Eliashberg equations 



is the self-energy matrix with its 



in V ' m.n=i 

D™ n (q,in) lm G rl (k F + q,iu; + iQ) ln + [ 
D% v (q, ity-f^G^ (k F + q,iw + iCl)j u T3, 



E ( — ^ tr l7mG^(fc + g, + ifi)7 n G„(fc, iw)]^ 

771,71 — 4 

^-itr[7 M T 3 G^(fc + g, iu + Olj^TaG^k, iw)]^ , (43) 
where this is basically the same as the previous one Eq. 



(17), but the holon propagator is replaced with a ma- 
trix including its superconducting part. If we rewrite the 
above equations with each component of the holon prop- 
agator, we obtain 

E(fc, iuj) 

4 

= m l T E E "t nG ( k + q,iu + i$l)D v nn (q, ifi) 7 „ 

iQ,q 71—3 

+T E E 7iG(fc + «,»w + in)X>&(?,in)7j (44) 

in : q i,j—x : y 

and 

A(k,iu) 

4 

= m l T E E T" F ( fc + «> iW + ift ) D »n(«» i0 )7« 
71—3 

- T E E + iw + tfWj (g, i0)7i, (45) 

in,q i,j—x,y 

where the normal and abnormal holon propagators are 
given by 



F(k,iw) 



g 1 (k,iuj) + H(k,iu>) 



[ 5 -!(fc, iw) + £(fc, jw)] 2 - I A(fc, iw)| 2 ' 
A(fc,zw) 

~ kr 1 ^, «w) + E(jfc, m;)] 2 - I A(fc, ito)\ 2 



(46) 



As shown in Eqs. (45) and (46), gauge fluctuations 
cause attractive interactions between holons with differ- 
ent isospins, guaranteed by the the presence of the T3 
matrix, while valance bond fluctuations result in repul- 
sive ones. The presence of repulsive interactions opens 
the possibility of d-wave pairing of doped holes. 

As discussed above, there is another pairing channel 
between holons with same isospins but on different sub- 
lattices, induced by valance bond fluctuations. Com- 
parison of these two superconducting phases by solving 
each Eliashberg equations remains as an important fu- 
ture work. What we would like to mention is that be- 
cause both gauge and valance bond fluctuations give rise 
to competing interactions, the present anomalous metal- 
lic phase is expected to be stable against superconduc- 
tivity in some parameter regions of the effective theory 
Eq. (10). Even if such a metallic phase becomes unstable 
against superconductivity at zero temperature, at finite 
temperatures the present non-Fermi liquid physics would 
survive. 

There is another instability channel associated with 
particle-hole pairing. Different from the particle-particle 
channel, gauge fluctuations cause attractive interactions 
for the particle-hole channel when holons have the same 
isospin. As pointed out by Altshuler et al., 3l| the 2kp 
vertex with the Fermi momentum kp is power-law diverg- 
ing with its exponent 1 /N approximately, where N is the 
flavor number of fermions. They showed that divergence 
of particle-hole susceptibility depends on the exponent 
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of the diverging vertex. In other words, the particle-hole 
pairing instability arises when the divergence of the 2kp 
interaction vertex is sufficiently strong. In this respect 
the present effective theory with the large N approxima- 
tion is expected to be free from such particle-hole insta- 
bilities. Remember that the large N approximation is 
consistent with the self-consistent Eliashberg framework. 
Interestingly, valance bond fluctuations are against such 
particle-hole pairings as they compete with gauge fluctu- 
ations for superconductivity. More quantitative analysis 
for the particle-hole channel is required when the flavor 
number of fermions is not large. However, this is be- 
yond the scope of the present paper since it is not clear 
whether even the Eliashberg approximation is stable or 
not in this case. 

In the recent publication [32j the presence of an inter- 
esting inhomogeneous superconducting state was demon- 
strated in the two-leg ladder system based on the renor- 
malized mean-field theory, where a heavy numerical anal- 
ysis was performed for solving self-consistent mean-field 
equations of order parameters in real space. From the 
present approach it is difficult to see the emergence of 
such complicated inhomogeneous order parameter pat- 
terns because our effective field theory is based on the 
uniform phase of order parameters. 

Applying the present theoretical framework to the two- 
leg ladder system, one would see that gauge fluctuations 
enhance the 2kp particle-hole vertex to cause power-law 
divergence, but the particle-hole susceptibility does not 
diverge at least in the large N approximation, as dis- 
cussed above. However, charge and spin density waves 
can certainly occur when the fermion flavor number is 
small. In this respect, when the coupling strength be- 
tween doped holes and valance bond fluctuations is large 
but the fermion flavor number is not large, an inhomoge- 
neous superconducting phase is expected to appear. Un- 
fortunately, fully self-consistent analysis including both 
particle-particle and particle-hole instabilities is clearly 
beyond the scope of the present paper. 



VII. SUMMARY 

Fermionizing the charge sector and bosonizing the spin 
part in the SU(2) slave-boson theory, we have derived 
an effective field theory for dynamics of doped holes 
in the antiferromagnetically correlated spin background, 
where spin fluctuations are described by the SO (5) 
WZW theory while charge dynamics is expressed by 
non-relativistic QED3 around four Dirac nodes. In par- 
ticular, hole dynamics affects deconfinement of bosonic 
spinons in the SO (5) WZW theory through the coupling 
term between valance bond (monopolc) fluctuations and 
fcrmionic holons. Such interactions give rise to z = 3 
criticality for monopole dynamics, prohibiting their pro- 
liferation in the presence of the WZW term. As a result, 
holon fluctuations turn out to help spin fractionalization 
near the quantum critical point. 



We have investigated thermodynamics and transport 
in the Eliashberg framework for our effective field the- 
ory. We pointed out that the Eliashberg framework is 
the lowest order self-consistent approximation well con- 
trolled in our effective field theory, where Migdal theorem 
works well owing to z — 3 criticality, and the large N v 
limit with the number of Dirac nodes N v is naturally al- 
lowed. We find that spin fluctuations are described by 
z = 1 for antiferromagnetic fluctuations and z — 3 for 
valance bond excitations, giving rise to three regimes, 
where superspin fluctuations are gapped at low temper- 
atures, only valance bond excitations are critical at in- 
termediate temperatures, and superspin fluctuations are 
critical at high temperatures. Both valance bond and 
gauge fluctuations are described by z = 3 critical theory, 
and we find non- Fermi liquid physics for thermodynamics 
and electrical transport near the quantum critical point, 
consistent with z = 3 scaling. In addition, even in the 
quantum disordered phase such non- Fermi liquid physics 
is preserved owing to critical gauge fluctuations. 

To further justify our effective field theory, we have ap- 
plied it to one dimension, physically well known. In one 
dimension spin fluctuations are described by the SO (4) 
WZW theory while charge excitations are represented 
by QED2. We have taken the monopole- holon coupling 
term into account in the abelian bosonization framework. 
We demonstrated that holon dynamics results in con- 
fining interactions between monopole excitations. Thus, 
we conclude that charge fluctuations help spinon decon- 
finement in both one and two dimensional cases while 
charge dynamics suppresses monopole fluctuations more 
strongly in one dimension. 

We have discussed stability of the non-Fermi liquid 
metallic phase against superconductivity and density 
waves. An interesting observation is that two kinds of 
holon pairing channels exist due to the presence of the 
isospin quantum number in the SU(2) slave-boson de- 
scription. In the different isospin channel attractive pair- 
ing interactions are caused by gauge fluctuations while 
in the same isospin channel such interactions arise from 
valance bond fluctuations. An important thing is that 
such gauge and valance bond interactions compete with 
each other. As a result, superconductivity is expected 
to appear in a limited parameter range. The presence 
of repulsive interactions allows the possibility of d-wave 
pairing of doped holes. For the particle-hole channel, 
we argued that as far as the fermion flavor number is 
sufficiently large, consistent with the Eliashberg approx- 
imation, the homogeneous metallic phase can be stable 
against charge and spin density waves. 

We would like to emphasize that the SU(2) structure 
is important for our treatment. If we start from the U(l) 
slave-boson representation, the fermionization procedure 
is not performed naturally since we have nonzero net flux 
owing to the presence of finite density of holons. Further- 
more, the SO (5) WZW theory for the spin sector does not 
appear owing to the contribution from a chemical poten- 
tial term in the U(l) slave-boson framework. 
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Our effective field theory exhibits direct interactions 
between monopoles and holes, associated with dcconfine- 
ment of bosonic spinons away from half filling. The im- 
portant issue how doped holes affect spinon deconfmc- 
ment deserves to be studied more carefully. 
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